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ABSTRACT 



We present Submillimeter Array observations of H2CO and N2H+ emission in 
the disks around the T Tauri star TW Hya and the Herbig Ae star HD 163296 
at 2"-6" resolution and discuss the distribution of these species with respect 
to CO freeze-out. The H2CO and N2H+ emission toward HD 163296 does not 
peak at the continuum emission center that marks the stellar position but is 
instead significantly offset. Using a previously developed model for the physical 
structure of this disk, we show that the H2CO observations are reproduced if 
H2CO is present predominantly in the cold outer disk regions. A model where 
H2CO is present only beyond the CO snow line (estimated at a radius of 160 AU) 
matches the observations well. We also show that the average H2CO excitation 
temperature, calculated from two transitions of H2CO observed in these two 
disks and a larger sample of disks around T Tauri stars in the DISCS (the Disk 
Imaging Survey of Chemistry with SMA) program, is consistent with the CO 
freeze-out temperature of ~20 K. In addition, we show that N2H+ and H2CO 
line fluxes in disks are strongly correlated, indicative of co-formation of these 
species across the sample. Taken together, these results imply that H2CO and 
N2H+ are generally present in disks only at low temperatures where CO depletes 
onto grains, consistent with fast destruction of N2H+ by gas-phase CO, and in 
situ formation of H2CO through hydrogenation of CO ice. In this scenario H2CO, 
CH3OH and N2H+ emission in disks should appear as rings with the inner edge 
at the CO midplane snow line. This prediction can be tested directly using 
observations from ALMA with higher resolution and better sensitivity. 

Subject headings: protoplanetary disks; astrochemistry; stars: formation; ISM: 
molecules; techniques: high angular resolution; radio lines: ISM 



- 3 - 



Introduction 



Planetary systems are assembled from dust and gas in the disks surrounding pre-main 
sequence stars. The nature of the formed planets are intimately linked to the structure, 
composition and evolution of the parent circumstellar disk. Molecular emission lines serve 
as probes of disk characteristics, such as density, temperature and ionization fraction, that 
are not accessible by other observations. For example, N2H+ along with the deuterated 
ions like H^ D+ and DCO''' are b elieved to trace the ionization fraction near the midplane 



of the disks ( lOberg et al. 



201ld ). Molecular distributions in disks are also important to 



characterize because of their connection with the composition of forming planetesimals. 
This especially true of organic molecules, of which H2CO is an important representative. 
While most molecules are expected to be reprocessed in larger planetary bodies, the 
disk-chemical composition m ay survive quite intact in icy planetesimals, including comets 



( iMumma fc Charnley 



20111 ). Such planetesimals may have seeded the Earth with water 
and organics, connecting disk chemistry with the origins of life. Predicting the organic 
composition of these planetesimals depend on our understanding of the distribution of the 
organic composition of grains in the disks. 

The overall disk chemical structure is set by a combination of photoch emistry at the 



disk surface and sequential freeze-out of molecules in the disk interior (e.g. 



Aikawa et al. 



2OO2I ). The distributions of N2H"'" and H2CO and their relationship to CO present important 
test cases of the chemical models. The CO molecule is one of the last to freeze out and 
predicted to deplete quickly from the gas-phase at T< 20 — 25 K for typical disk mid-plane 
densities. Abundant N2H+ is expected only where CO is depleted because N2H"'" forms from 
prot onation of No, whi ch remains in the gas-phase at temperatures a fe w degrees lower tha n 



CO flOberg et al. 



2OO5I ). and is destroyed mainly by reactions with CO (IBergin et al.ll2002l ). 



In dense cores in star forming regions, the abundance of CO is observed to show a strong 
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anti-correlation with N2H"'" (e.g. ICaselli et al.lll999l : iBergin et al.ll2002l : |j0rgenseru |2004| ) . In 
disk models, this effect is manifested as a jump in the N^H + column density at the CO 



2012h . The snow 



"snow line" , one order of magnitude in a recent calculation (iWalsh et al.l 
line is here defined as the disk radius where the midplane dust temperature is cold enough 
for volatiles to condense into ice grains. 

The H2CO chemistry is more complicated than N2H"'" because H2CO can form through 
several different pathways, both in the gas-phase and on grain-surfaces. Grain-surface 
formation of H^CO should depend directly on CO freeze-out; con strained by theory 



(iTielens fc 



iagen 



Fuchs et al. 



1982 



Cuppen et al 



Watanabe et al. 



2003 



20091 ) and experiments (e.g. 
20091 ). H2CO (and CH3OH) form readily from CO ice hydrogenation. If this 
is the dominant for mation pathw'ay of H^CO in disks, and if ices are partially desorbed 



non-thermally (e.g. 



Garrod et al. 



20071 : 



Oberg et al 



2009b 



a|), then H2CO gas should 



coincide with N2H+ exterior to the CO snow line. 

To date, emission from millimeter wavelength H2CO and NoH"*" lines has been detecte d 



Oi et al. 


2003 


201J, 


2OI1J). 



Thi et al. 



2004; 



disks, respec 


ively 


Dutrev et al. 


2007 



Dutrev et al. 



1997 



Aikawa et al 



Henning fc Semenov 



2008 



2003 



Oberg et al. 



Most detections are toward T Tauri stars with massive disks, and the 



201 Ibl ). Based 



detection fraction toward more luminous Herbig Ae stars is low ( lOberg et al.l 
on these observations, H2CO and N2H+ are mainly abundant in disks with large reservoirs 
of cold dust and gas, where CO freeze-out is expected to occur. A direct connection between 
CO freeze-out and N2H+ and H2CO in disks has yet to be observationally established, 
however. 

In this paper we present Submillimeter Array (SMA) observations of H2CO and N2H+ 
toward the disks around HD 163296 and TW Hya, and we use these new observations along 
with H2CO and N2H"^ observations from DISCS (Disk Imaging Survey of Chemistry with 
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SMA) to constrain the H2CO and N2H+ distributions. The new data and their cahbration 
are described in fJ2J In §21 we present the H2CO and N2H"'" images and spectra toward 
HD 163296 and TW Hya, models of the H2CO distribution toward HD 163296, H2CO 
excitation temperature calculations, and examine the relationship between H2CO and N2H"'" 
emission across the sample of disks. In ^ we discuss the implications of these results, 
summarize the mounting evidence for CO-ice regulated chemistry and make predictions 
for future observations of H2CO and N2H"'" emission from disks with better sensitivity and 
resolution. 

2. Observations 

The observations of HD 163296 (R.A. = 17^56'"21.279% decl. = -2r57'22'.'38; 
J2000.0) were made between 2008 and 2012, and of TW Hydrae (R.A. = ll'^01'^51.875^ 
decl. = -34°42'17'.'155; J2000.0) between 2008 and 2012, using the eight-antenna 
Submillimeter Array (SMA) located atop Mauna Kea, Hawaii. Table [1] provides a summary 
of the observational parameters and results. For the 2007 and 2008 observations, the SMA 
receivers operated in a double-sideband mode with an intermediate frequency (IF) band of 
4-6 GHz from the local oscillator frequency, sent over fiber optic transmission lines to 24 
overlapping "chunks" of the digital correlator. The 2012 observations were made after an 
upgrade that enabled a second IF band of 6-8 GHz, effectively doubling the bandwidth. 

The SMA observations of HD 163296 were carried out in the compact-north (COM-N), 
compact (COM) and subcompact (SUB) array configurations. The 2007 observations 
included the DCO+ 3-2 line at 216.1126 GHz and the H2CO 3i,2 - 2i,i line at 225.698 GHz. 
The 2012 observations included the N2H+ 3-2 at 279.512 GHz and the H2CO 4i,4 - 3i,3 
line at 281.527 GHz. The observing loops used J1733-130 as the main gain calibrator and 
observed J1744-312 every other cycle to check the phase calibration. Flux calibration was 
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done using observations of Titan and Uranus. The derived fluxes of J1733-130 were 1.19 
Jy (2007 Mar 20), 1.25 Jy (2012 Jun 10), 1.40 Jy (2012 Aug 12 and 14). The bandpass 
response was cahbrated using observations of 3C279, Uranus and J1924-292. 

The SMA observations of TW Hya were carried out in the compact (COM) and 
subcompact (SUB) array configurations. The 2008 observations included the H2CO 
5i,5 - 4i,4 line at 351.769 GHz. The 2012 Jan 13 SUB observation included the H2CO 
4i,4 — 3i,3 and N2H+ 3-2 lines, like HD 163296, but unfortunately the chunk containing 
N2H+ was corrupted and unusable. The 2012 Jun 04 COM observation using a similar 
setting successfully included the N2H"'" line. The observing loops used J1037-295 as the gain 
calibrator. Flux calibration was done using observations of Titan and Callisto. The derived 
fluxes of J1037-295 were 0.73 Jy (2008 Feb 23), 0.73 Jy (2012 Jan 13) and 0.82 Jy (2012 
Jun 4). The bandpass response was calibrated using observations of 3C279 and 3C273. 

Routine calibration tasks were performed using the MIR software package Q, and 
imaging and deconvolution were accomplished in the MIRIAD software package. 



3. Results 

In this section we present detections of II2CO and N2H"'" emission lines toward 
HD 163296 and TW Hya, display their respective distributions ( §3.ip . and compare the 
higher quality observations of H2CO in HD 163296 with models ( §3.2p . We then combine 
the new data with previously reported H2CO and N2H'^ detections in disks to examine 
trends with respect to H2CO excitation temperature ( §3.3p and with each other ( §3.4p . 



^http:/ /www. cfa.harvard.edu/~cqi/mircook. html 
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3.1. HaCO and N2H+ towards HD 163296 and TW Hya 



Figure [T] shows images of the spectrally integrated emission toward TW Hya and 
HD 163296 at the rest frequencies of two H2CO lines and the N2H+ J = 3 - 2 line. H2CO 
4i,4 — 3i^3 and 5i^5 — 4i 4 and N2H+ are detected toward TW Hya, and H2CO 3i^2 — '^1,1 and 
H2CO 4i,4 ~ 3i^3 and N2H^ are detected toward HD 163296. The emission toward TW Hya 
appears to be centrally peaked at the size scale of the beams (FWHM > 2"). 

Toward HD 163296, however, neither the H2CO lines nor the N2H"'" line emission peaks 
at the location of the continuum peak that marks the stellar position, but instead show 
significant offsets. This is most readily apparent for the H2CO 3i^2 — 2i 1 line that was 
observed with a slightly smaller and more advantageously rotated beam, where the emission 
appears ring-li ke. This is the second repo rted observation of a ring-like H2CO distribution 



after DM Tau (IHenning fc Semenov 



20081 ). Interpreting the H2CO emission toward DM Tau 



is however complicated by a large central cavity in dust emission (lAndrews et al. 



201 ih . 



Figure [2] shows the spatially integrated spectra. The line shapes and central velocities 
agree with what has been previously observed for other molecular lines toward these disks. 
The line fluxes are listed in Tabled], and the values are comparable to detections of these 



lines toward other large protoplanetary disks (lOberg et al. 



2010 



2011br i. 



3.2. HD 163296 H2CO Model Results 

Figure [T] demonstrates that H2CO emission toward HD 163296 is spatially resolved 
and thus contains information on the radial distribution of H2CO in the disk. The relative 
excitation of the two H2CO transitions should probe primarily the vertical distribution 
and thus provide complementary constraints. We explore the H2CO distribution based on 
a previously developed accretion disk model with a well-defined temperature and density 



- 8 - 



structure, constrained by the HD 163296 broadband spectral energy distribution, spatially 
resolved millim eter dust continuum, and multiple CO and CO iso topologue line observations 



( lOi et al 



Qi et al. 



( 120081 ) for constraining the H2CO 



20111 ). We adopt the same methods as 
distribution, here fitting models that assume a radial power-law ( §3.2. ip and a simple ring 
with inner boundary at the CO "snow line" ( §3.2.2^ . 



3.2.1. Power-law Model 



For a first-order analysis of the distribution of H2CO, we model the radial variation in 
the column density as a power law Nioo x (r/100)^ between an inner radius Rj„ and outer 
radius Rout, where Nioo is the column density at 100 AU in cm~^, r is the distance from the 
star in AU, and p is the power-law index. 

For the vertical distribution, we assume that H2CO is present with a constant 
abundance in a layer with boun daries toward the midplane and toward the surface of the 



disk (similar to 



Oi et al. 



Aikawa fc Nomura 



( I2OO8I )). This assumption is motivated by chemical models (e.g. 
2006f ) that predict a three-layered structure where most molecules are 
photodissociated in the surface layer, frozen out in the midplane, and have an abundance 
that peaks at intermediate disk heights. The surface (cr^) and midplane (o"m) boundaries 
are presented in terms of S21 = Sj://(1.59 x lO^^cm^^), where Hh is the hydrogen column 
density measured from the disk surface. This simple model approximates the vertical 
location where H2CO is most abundant. The excitation of multiple transitions can constrain 
both as and am, but in this case of v e ry mo dest signal-to-noise, we fix as to 0.79, the surface 



boundary found for CO by 



Qi et al. 



teOllI ). and we fit a^n for the midplane boundary and 



the power-law parameters (Nioo, P, R-m and R, 



out ) 



Using the structure model, we compute a grid of synthetic H2CO visibility datasets 
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over a range of Rout, R-in, P, o'm and Nioo values and compare with the observations. The 
best-fit model is obtained by minimizing x^, the weighted difference between the real 
and imaginary part of the complex visibility measured in the {u, f )-plane sampled by the 
SMA observations of both H2CQ transitions. We u se the two-dimensional Monte Carlo 



model RATRAN (iHogerheijde &: van der Tak 



20001 ) to calculate the radiative transfer and 



molecular excitation. The coUisional rates are taken from the Leiden Atomic and Molecular 



Database (ISchoier et al. 



20051 ). 



Table [2] lists the best-fit parameters of the model. The power-law index of 2 implies an 
H2CO column density that strongly increases with radius. Figure |3] shows the best-fit radial 
distribution of the H2CO column density. Figures H] and [5] present comparisons between the 
observed channel maps and the best-fit model. The model reproduces the main features of 
the observations remarkably well, in particular the fiux ratio between the inner and outer 
channels, and the lack of emission at the location of the continuum peak. 

Figure [6] shows the surfaces for the Rj„ and Rout versus the power law index p, 
which enables us to quantify the uncertainties associated with the inner and outer region 
sizes and the power-law index. We find that p is constrained between 0.5-3.0 (within la) 
while Kin is constrained to be <200 AU. The value does not change significantly for 
inner radii <90 AU, as expected for the 2" beam size of the observations. The outer radius 
is better constrained, since the emission is very sensitive to the value of Rout with a positive 
power-law index. Figure [7] shows the H2CO 3-2 line spectrum compared with the spectra 
derived from models with different radial column densities power-law indices. The spectra 
suggest a lack of high velocity line wings associated with emission originating in the inner 
regions of the disk, consistent with the results of the analysis. 
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3.2.2. Ring Model 

The positive power-law index found for the H2CO radial distribution implies that 
H2CO is present mainly in the outer disk. This is expected if H2CO forms in situ from 
CO ice hydrogenation and is th erefore pr e sent m ainly beyond the CO snow-line, previously 



determined to be at 160 AU by 



Oi et al. 



(120111 ). Guided by this astrochemical ansatz, we 



have tried a second "ring" model where the H2CO gas is only present where CO has frozen 
out. The vertical su rface boundary is then defined by the CO freeze-out temperature of 



19 K flOi et al. 



2OIII ). while the midplane boundary can be constrained by the excitation of 
multiple H2CO transitions, as in the power-law model. Within this layer, the abundance of 
H2CO is assumed to follow that of H nuclei with a constant fractional abundance, which is 
also a parameter fit to the data. 

The best-fit abundance is 5.5 x 10^^^ and the midplane boundary a-m is consistent with 
what we find in the power-law model (Table [2]). Figure [3] shows that the vertical surface 
boundary at 19 K effectively results in a ring-like radial structure, where the inner edge of 
the ring is at CO snow line. Figure |3] also shows that the best-fit power-law and ring models 
result in similar H2CO column densities beyond the CO snow line. The profile of the ring 
model is considerably fiatter than the power-law model and even drops outside of 300 AU 
exponentially to the edge of CO emission. Figures HHS] show that the power-law model and 
the ring model channel maps display some subtle differences. But both provide good fits to 
the data within the noise of the SMA observations. The same model also provides a good 
match to the N2H"'" fiux, but the combination of low signal-to-noise and observations of just 
one transition preclude any independent modeling of the N2H"'" distribution. 



3.3. H2CO Excitation Temperatures 



HD 163296 is very favorable for H2CO and N2H+ imaging since its relatively high 
luminosity and massive disk puts the CO snow line at a large angular distance compared 
to other disks, enabling us to resolve the H2CO and N2H"'" emission. Without such spatial 
information, however, we can still obtain a constraint on where H2CO emission originates 
in disks based on the average H2CO excitation temperatures. As a gross approximation, 
H2CO that coexists with CO ice should be cold, i.e. present at an excitation temperature 
comparable to the CO freeze-out temperature. To test the viability of using excitation 
temperatures to constrain the H2CO distribution we extracted spectra from three of the 
HD 163296 simulations presented in §3.2, selecting models with the best outer radius, 
no inner hole and power- law indices of —2, 0, and 2. These distributions approximately 
correspond to a H2CO abundance that follows the H2 column, that keeps constant with 
radius and that increases steeply, forming a ring. Assuming LTE, that both H2CO 
transitions trace the same underlying populations, and a single rotational excitation 
temperature, T^^t, we calculate 

™* " In {{v,Si4jT^dv)l{voSi4!Trdv)) ' 

where Eq and Ei are the upper energy levels for the low and high H2CO transitions used in 
the calculation (H2CO 3o,3 — 2o,2 and 4i^4 — 3i^3 for most disks), v and Sji^ the corresponding 
line frequencies and temperature independent transition strengths and dipole moments, 
J Tdv the integrated line intensity, which is calculated from the integrated fluxes based on 
F/T — 13. 6A^/ (a X 6), where F is the flux in Jy, T the intensity in K, A the line wavelength 
in millimeters, and a and h the emission diameters in Because of both vertical and radial 
temperature gradients in the disk, the size dependence of the emission regions from two 
transitions is complicated but the emission area is not expected to be very different. For 
simplicity we assume the extent of the emission is the same for both transitions on account 
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of the model dependent effects of the temperature gradients. All of the line parameters 
were gathered from Splatalogue (a transition-resol ved compilation of several spectroscopic 



databases), with the data originating from CDMS ( iMiiller et al. 



20051 ) ■ LTE is a reasonable 



approximation if H2CO is mainly present at high densities. Th e critical densities for the 



observed H2CO transitions vary between 10^ and 5x10^ cm ^ ( iTroscomptet al. 



20091) 



dependent on assumed kinetic tem peratures. At radii <300 AU, all gas colder than 25 K is 



at densities higher than 10^ cm ^ (iQi et al 



201 ll ). justifying this assumption. 



Using the simulated spectra we derive an excitation temperature of 27 K for the model 
with p = —2 and excitation temperatures close to or below 20 K for the other two models. 
Excitation temperatures thus provide some constraints on the H2CO distribution, but a 
temperature close to the expected CO freeze-out temperature only implies an increasing 
abundance with radius; it cannot be used to assess how steeply the abundance increases. 

By combining the new H2CQ detection s toward TW Hya and HD 163296 with H2CO 



data from DISCS (jOberg et al. 



2010 



2011bl ). we have a sample of 10 disks with two H2CO 
line detections or one H2CO line detection and one upper limit (Table [3]). These data are 
sufficient to calculate excitation temperatures, albeit with substantial uncertainties. Figure 
[8] shows the calculated excitation temperatures for 9 of the 10 disks; the chemically peculiar 
Herbig Ae star HD 142527 is not included because its excitation temperature of 250 K is 
probably not due to thermal excitation. 

The H2CO excitation temperature, listed in Table [3], is consistent with, or lower 
than, the CO freeze-out temperature of ~20 K for all of these disks. The average H2CO 
temperature in the sample (excluding HD 142527) is 18 ± 6 K. We note that the relatively 
high excitation temperature toward HD 163296 is most likely due to the fact that some 
H2CO 3-2 emission is resolved out by the SMA observations, as the compact-north antenna 
configuration has few short baselines. 
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3.4. H2CO/N2H+ Correlations 

We examine tlie disk sources to test if N2H"*" and H2CO emission are correlated 
across the sample, as would be expected if 1) the two molecules form under the same 
physio-chemical conditions, i.e. only in the regions where CO has frozen out, 2) the line 
emission trace the total N2H+ and H2CO column well, and 3) midplane ionization levels 
and CO hydrogenation efficiencies do not differ 'too much' across the sample (i.e. the size 
of the CO freeze-out region is the most important regulator of N2H+ and H2CO column 
across the sample). It should be noted that there are other scenarios that could produce 
a correlation as well, and that the correlation analysis below should only be considered 
as a constraint on the H2CO distribution in combination with the results in the previous 
sections. In particular a constant H2CO/N2H+ across a disk sample would be expected if 
the relative fractions of chemically characteristic disk regions is always similar in disks. To 
conclusively test this requires a larger sample than currently available, but the fact that we 
did not find that H2CO emission correlates with any other molecular emission than N2H^ 
already challenges this scenario. 

Where possible, we base the comparison on the H2CO 3o,3 — 2o,2 line that has 
-Eup = 21 K, similar to N2H+ J = 3 — 2 {E^p = 27 K), to minimize variations in fiuxes 
caused by the different detailed temperature structures in different disks. Excluding HD 
142527, this line has been observed toward 6/9 of the sample. For the remaining 3/9 
disks, we calculate the expected H2CO 3o,3 — 2o,2 line flux based on the H2CO excitation 
temperatures and fluxes of other H2CO lines toward each source. We then normalize 
the flux of each H2CO 3o,3 — 2o,2 line to the (Taurus) distance of 140 pc, and we further 
normalize to a disk mass of 0.01 Mq to account for the fact that more nearby and more 
massive disks tend to have overall stronger line emission. Figure |9] shows that there is a 
strong correlation between the normalized H2CO and N2H+ fluxes in the disk sample; the 
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rank correlation is statistically significant at the 95% level. As expected, Figure [TO] shows 
that this implies a nearly constant N2H+ 3-2 / H2CO 89,3 — 2o,2 flux ratio across the sample. 



4. Discussion 



( 120081 ) has been to first constrain 



Our modeling strategy in this study and in iQi et al. 
the overall structure of molecular emission in disks using a parametric model with a 
minimum of free parameters, i.e. to determine whether the radial column density profile 
of a species decreases , increases or is fiat as a function of disk radius. As demonstrated 



m 



Oberg et al. 



(120121 ). the slope of the radial column density profile already can place 
significant constraints on the formation pathway of a molecule. Here we find that H2CO 
toward HD 163296 belongs to the family of molecules that display an increasing column 
density with radius. This first-order constraint on the H2CO distribution motivated us to 
consider H2CO formation pathways that would result in an increase of H2CO with disk 
radius. H2CO formation through CO-ice hydrogenation results in a simple prediction that 
H2CO should be present in a ring, with the inner edge at the CO snow line. We therefore 
set up a second model, based on this prediction, to test if the observations are consistent 
with this hypothesis for H2CO formation. We propose that this combination of backward 
and forward modeling both provides a fair view of the constraints obtained by fitting the 
data, and challenges our basic understanding of disk chemistry. 



4.1. H2CO Formation 



H2CO can form through multiple chemical pathways. We have shown that the H2CO 
distribution towards HD 163296 and the sample statistics are consistent with formation 
through in situ CO ice hydrogenation. Here we consider the effects of additional pathways 
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for H2CO formation, in particular (1) in situ gas phase formation, and (2) formation in the 
pre- and proto-stellar phases, followed by incorporation into the disk. 

H2CO can form in the gas-phase through ion-neutral r eactions involving, e.g. C Hj" or 



through neutral-neutral reactions between CH3 and O (e.g. 



Aikawa fc Herbst 



19991). The 



neutral-neutral formation pathway is expected to result in a radially flat column density 



structure for a typical T Tau ri disk (lAikawa fc Herbst 



at temperatures of 20-40 K ( lAikawa et al 



I999I ). with most emission originating 



2OO3I ). It is not clear from existing disk models 



whether the structure will look sub stantially different if th e neutral-ion reactions dominate 



In a protostellar chemistry model (IBergin &: Langer 



19971 ). CH3 disappears when CO 



depletes (E. Bergin, private communication). This suggests that H2CO forming through 
this pathway should be anti-correlated with CO-freeze-out. Neither of these gas-phase 
pathways thus predicts excess H2CO column densities in the outer disk, or at low (T<20 K) 
temperatures. The observed low excitation temperature could on its own be explained by 
efficient turbulent mixing of H2CO formed in the gas-phase and then cooled down in the 
midplane regi ons, similarly to what has been proposed to explain cold CO gas in disks 



( lAikawa 



20071 ). Turbulent mixing would not, however, explain the observed deficiency of 
H2CO towards the inner disk in HD 163296. Gas-phase formation of H2CO then seems an 
unlikely dominant source of H2CO in disks in light of the new observations, but the case 
is unlikely to be conclusively settled until the exclusive grain-surface product CH3OH is 
observed to display a similar distribution. 

H2CO in disks could be a product of protostellar or molecular cloud chemistry, as 
H2CO is commonly observed in pre- and protostellar sources, and this molecular content 
may be preserved, at least in part, through the process of disk formation, 
has modeled this scenario, starting with a H2CO ice abundance of lO^^rin inherited from 
the cold cloud. The model also includes H2CO formation through gas-phase processes 



Willacvl (1200 
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in the disk and results are presented with and without photodesorption. In the model 
without photodesorption, H2CO follows CO in the inner disk and has an additional outer 
disk component with an abundance that decreases with radius beyond the CO snow line. 
When photodesorption is included, the H2CO column is flat across the disk, corresponding 
to a power-law index of 0. Neither predicted abundance pattern is consistent with the 
new observations. I n addition, HoCO i s common towards pr otostellar sources of a range 



of luminosities (e.g. 



Schoier et al 



2004 : 



Bisschop et al 



20071 ). while it is pre- dominant ly 



detected towards disks around the low luminosity T Tauri stars. 

In short, II2CO formation through in situ CO-ice hydrogenation is not only consistent 
with the observations, but it is the only pathway proposed (so far) that naturally explains 
the observations. To conclusively demonstrate a CO-ice hydrogenation origin would, 
however, require the detection of co-spatial emission of CH3OH; CH3OII has no known 
efficient gas-phase formation path way and is predicte d to form together with H2CO 



whenever CO ice is hydrogenated ( ICuppen et al 



20091). 



4.2. Locating the CO "Snow Line" 



The "snow line" is typically used to denote the midplane disk radius at which the 
temperature is low enough for water to condense out on dust grains. Outside of the 
snow line, grain accretion will be faster because of larger and stickier grains, whi ch may 



substantially speed up the f ormation of planetesimals and eventually planets (e.g. 



1981 



Ida fc Lin 



2004 . 



2008 



Ciesla fc Cuzzi 



2006 



Kretke fc Lin 



Hayashi 



20071 ). In the Solar System, 



the diyiding line between rocky planets and gas giants coincide with the H2O snow line 



([Lewis 



19741 ). CO is another abundant volatile in disks and its snow line may boost plane t 



formation in the outer disk by providing extra solid masses (IDodson-Robinson et al. 



and inducing planet traps flMasset et al 



2006 



Hasegawa &: Pudritz 



20091 ) 



2OI2I ). and could affect 
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the elemental make-up of the forming gas-giants (lOberg et al.ll2011al ). Because of its 
high volatility (Tfreezc-out ~ 20 K), the CO snow line is expected at disk radii of lOs-lOOs 
of AU. This should make it a far more accessible target than the H2O snow line for 
millimeter interferometry studies aimed at examining the general effects of snow lines on 
disk structures. 

Localizing the CO snow line directly from millimeter CO data is challenging, however. 



from the central 


star, £ 


lAikawa & Herbst 


1999) 



one set by radiative heating at the disk surface (e.g. lAikawa fc HerbstI 
in a CO condensation front that is located at different radii at different disk heights (Fig. 
[TT]) . and also that some CO is present at all radii in the upper disk layers. This fact, 
together with a complex radiative transfer (most CO lines are expected to be optically 
thick in the disk center and optically thin in the outer parts of the disk), means that 
the location of a CO snow line in a disk cannot be inferred from simply inspecting a CO 
disk image. The best constraint that exists to date on a CO snow line radius is toward 
the Herbig Ae star HD 163296 based on the analysis of multi-transition, multi-isotope , 



20111). 



spatially resolved CO line data on a self-consistent physical disk model ( lOi et al. 
The temperature structure of the model has been constrained by optically thick multiple 
CO lines and detailed analysis of the optically thinner ^'^CO emission reveals a significant 
column density reduction at around 165 AU that cannot be explained by the overall disk 
column distribution as traced by the dust. This is interpreted as the result of CO freeze-out 
and the location as the CO snow line. Uncertainties in the temperature structure will 
mainly affect the determination of the CO freeze-out temperature at the location of the CO 
snow line, rather than the location itself , which is constrained between 135 and 175 AU 
in the disk of HD 163296. While fruitful, this is a time consuming approach that will be 
difficult to apply to larger samples of disks and will always include some degree of model 
dependency. 



- 18 - 



Another approach to constrain the CO snow hne location is to identify trace species 
that are only present where CO has begun to freeze out. Such molecules should display a 
ring-like structure with the inner edge corresponding to the midplane CO snow line. As we 
have described, N2H"'" and H2CO are good candidate probes of CO snow lines. In principle, 
these species can be used as powerful chemical imaging tools to constrain CO snow line 
locations in large samples of disks, rather than relying on complex analysis of the CO 
isotopologue observations. 

4-2.1. Simulated ALMA Observations 

It seems clear that both N2H"'" and H2CO are outer-disk species from the chemical 
perspective. To connect their formation to the onset of CO freeze-out conclusively requires 
a combination of higher sensitivity, and higher resolution imaging. This kind of imaging 
can be done with the newly available capabilities of the ALMA telescope in Chile, which is 
nearing completion of construction. 

To demonstrate the astrochemical predictions generated by our analysis of the SMA 
data, and the ease at which they can be tested with ALMA, we present a set of simulations 
of HD 163296 ALMA observations using the antenna configuration 5 (corresponding to 
0.3" resolution) in Figure [T2j The predicted II2CO and N2H+ rings are readily detected 
at this resolution, and the power-law and ring models are clearly distinguished. For both 
models, the diameter of the emission maximum can be estimated directly from the images 
within a fraction of the beam size, enabling us to determine if the II2CO and N2H+ emission 
truly trace the CO snow line. In addition, ALMA should have sufficient sensitivity to detect 
CH3OH if it is present with a similar abundance and distribution as H2CO, as expected 
if both of these species are formed through CO hydrogenation and then non-thermally 
desorbed. 
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5. Conclusions 

We have presented three observational results that support the idea that CO freeze-out 
regulates the H2CO and N2H"'" chemistry in disks: 

1. H2CO and N2H"'" emission towards HD 163296 appears offset from the continuum 
peak at a size scale consistent with the CO snow line at 160 AU. These observations 
are matched well by a ring model where H2CO is present only in the disk regions with 
CO freeze-out. 

2. The H2CO excitation temperature in a sample of 9 disks is typically below ~ 20 K, 
consistent with the bulk of H2CO emission originating in disk regions where CO is 
expected to freeze-out. 

3. The H2CO and N2H+ emission are correlated across the disk sample, consistent with 
the hypothesis of coexistence beyond the CO snow line. 

These results suggest that both N2H"'" and H2CO should be present in rings, with 
the inner edge at the CO snow hue. This may be used as a probe of CO snow hue 
locations across samples of disk and is also important for predicting the organic content of 
comets forming at different disk radii. In general, the radial and vertical distributions of 
molecules constitute strong probes of the basic chemistry used into astrochemical models, 
while molecular abundances and column densities are probably best used to test our 
understanding of the structure and history of individual objects. 

Facilities: SMA 

The SMA is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian 
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Fig. 1.— H2CO and N2H+ observations toward TW Hya and HD 163296 with the SMA. 
The emission toward TW Hya is centrally peaked, while there is a clear offset in the H2CO 
emission toward HD 163296. The first two contour levels are 3 and 5a, with la measured 
to be 0.10 (H2CO 4-3), 0.13 (H2CO 5-4), and 0.36 (N2H+) Jy km per beam toward 
TW Hya, and 0.06 (H2CO 3-2), 0.17 (H2CO 4-3), and 0.13 (N2H+) Jy km s^^ per beam 
toward HD 163296. 
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Fig. 2.— Spatially integrated spectra of H2CO and N2H+ toward TW Hya and H D 163296 



The r ed da shed lines mark Vlsr toward each source, based on CO observations (lOi et al. 



2006 



2OIII ). The double-peak structure typical for rotating disks is not resolved toward 



TW Hya with the applied spectral resolution because of its face-on orientation. 
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Fig. 3. — The hydrogen nuclei column density toward HD 163296 (iQi et al.l 
together with the H2CO column density from the best-fit power law and 'ring' models. In 
the ring model, the H2CO abundance is defined to be zero when the temperature is less than 
19 K, the CO freeze-out temperature, which leads to the sharp drop in the H2CO column 
density interior to the CO snow line at 160 AU. 
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Fig. 4. — Observed and simulated channel maps for H2CO 3i 2 — 2i 1 toward HD 163296, 
using the best-fit power-law (p=2) and ring model. Both models fit the data within the 
uncertainties. The first contour is 2a and each following contour step is la. The channel 
velocity in km s~^ is in the upper right corner of each panel and the synthesized beam is 
displayed in the lower left panel. 
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Fig. 5.— As Fig. Hbut for the H2CO 4i,4 - 81,3 line. 
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Fig. 7. — Simulated spectra for H2CO 3i,2 — 2i,i for models with radial column densities 
power-law indices p = 2(best-fit, red line), 0.5(within la noise level, green line) and — 2(blue 
line), overlaid with the HD 163296 spectra in gray shade. 
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Fig. 8. — The calculated H2CO excitation temperatures for all disks observed with the SMA 
that have at least one H2CO detection and one upper limit, except for the disk around HD 
142527, which has a very high excitation temperature of >100 K. All remaining disks have 
excitation temperatures H2CO consistent with or lower than the expected CO freeze-out 
temperature of ~20 K (red dashed horizontal line). 
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Fig. 9.— Correlation of N2H+ 3-2 (E^ = 27 K) and H2CO 3o,3 - 2o,2 (E^ = 21 K) emission 
normalized to disk mass (based on dust modeling) and source distance. The red symbols 
mark disks where the H2CO 3o,3 — 2o,2 flux has been calculated based on the flux from other 
H2CO transitions. 
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Fig. 10. — The N2H"'"/H2CO ratio, demonstrating that it is almost constant, as would be 
expected from the strong correlation between the normalized fluxes. 
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Fig. 11. — Illustration of the expected distribution of CO in the gas-phase and on grain- 
surfaces. Co freeze-out beyond the CO-snow line (R20K) and interior to the dashed contour 
is predicted to result in a large increase of N2H"'" and the onset of H2CO production from 
CO ice. The ice may then non-thermally desorb to produce gas-phase H2CO. 
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Fig. 12. — The predicted morphology of H2CO 3-2 emission toward HD 163296 for the 
power-law and ring model, when using ALMA antenna configuration 5 (corresponding to 
0.3" resolution) and Ih integration. The H2CO ring structure should be accompanied by 
similar CH3OH and N2H"'" rings if the observations reported in this paper are due the ex- 
clusive presence of H2CO in the outer disk, exterior to the CO snow line. For the CH3OH 
simulation we assumed the same column density and excitation temperature for CH3OH as 
has been observed for H2CO, and imaged the 241.767 GHz 5o,5 — 4o,4 line. White contours 
are [0.025,0.05,0.1] Jy km s^^ per beam. 



Tabic 1: Observational Parameters" 



Parameter 




HD 163296 






2007 Mar 20 (COM-N) 


2012 Jun 10 (COM) 


2012 Jun 10 (COM) 








2012 Aug 12 (SUB) 








2012 Aug 14 (SUB) 


Baseline (m) 


16.4-139.2 


16.4-77.0 


9.5-77.0 


Lines 


H2CO 3i,2-2i,i 


H2CO 4l,4-3l,3 


N2H+ 3-2 


Rest frequency (GHz) 


225.69778 


281.52693 


279.51170 


Beam Size (FWHM) 


3'.'4x2'.'2 


3^'9x2'.'5 


4'.'9x3'.'4 


P.A. 


54° 


15° 


14° 


Channel spacing (kms^^) 


0.54 


0.87 


0.60 


RMS Noise (Jybeam-i) 


0.063 


0.096 


0.072 


Integrated Flux (Jykms~^) 


0.89 


1.55 


0.53 


Continuum (GHz) 


221 


273 


273 


Beam Size (FWHM) 


3'.' 5x2'.' 2 


4yix2'.'6 


5'.'2x3y5 


RA. 


56° 


15° 


14° 


RMS Noise (mjybeam^^) 


3.4 


9.4 


3.8 


Flux Density (Jybeam^^) 


0.49 


0.96 


0.95 


Integrated Flux (Jy) 


0.61 


1.06 


1.07 


Parameter 




TW Hya 






2012 Jan 13(SUB) 


2008 Feb 23 (COM) 


2012 Jun 4 (COM) 


Baseline (m) 


9.5-45.2 


16.4-77.0 


16.4-77.0 


Lines 


H2CO 4l,4-3l,3 


H2CO 5l,5-4l,4 


N2H+ 3-2 


Rest frequency (GHz) 


281.52693 


351.76866 


279.51170 


Beam Size (FWHM) 


9'.' 1x5'.' 6 


3y3x2'.'0 


3'.'5x2'.'0 


P.A. 


24° 


1° 


10° 


Channel spacing (kms~^) 


0.43 


0.17 


0.11 


RMS Noise (Jybeam-i) 


0.16 


0.23 


0.55 


Integrated Flux (Jykms~^) 


1.22 


0.54 


2.2 


Continuum (GHz) 


273 


346 


273 


Beam Size (FWHM) 


9'.' 2x5'.' 7 


3^'5x2'.'0 


4'.'8x2{'4 


P.A. 


25° 


-2° 


10° 


RMS Noise (mjybeam^^) 


6.2 


12 


14 


Flux Density (Jybeam^^) 


0.90 


1.20 


0.70 


Integrated Flux (Jy) 


0.93 


1.52 


0.92 



"^AU quoted values assume natural weighting. 
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Table 2: Best-Fit Model Parameters^ 



Paranieter 


HD 163296 


Power-law Model 


Nioo (cm-2) 


1.7x10^^ 


P 


2 


Rout (AU) 


410 


Rin (AU) 


<200 


Surface boundary Us 


10-^-^^0.79 


Midplane boundary Um 


10^-5=31.6 


x' 


720101 


Ring Model 


Fractional abundance 


5.5x10-1^ 


Rout (AU) 


500 


Surface boundary 


T<19 K 


Midplane boundary cr^ 


101-5=31.6 




720101 



^Parameters in italics were fixed from CO modeling (cts for the power-law model and Rout 
for the ring model). 
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Table 3. Central star and disk, and H9CO data. 



Source 


RA 


DEC 


Spec. 


d 




Detected H2CO' 


Calc. 3o,3 - 2o,2 


Trot 








type 


[pc] 


[Mq] 




[Jy km s"-'] 


[K] 



TW Hya^' 


11 


01 


51.91 


-34 


42 


17.0 


K7 


56 


0.005 


4l,4 


-3i 


3, 5l,5 


- 4l,4 


~1.3j 


11(6) 


HD 163296'= 


17 


56 


21.29 


-21 


57 


21.9 


Al 


122 


0.089 


3l,2 


-2i 


1, 4i,4 


-3l,3 


~0.46j 


29(10) 


DM Tau<: 


04 


33 


48.73 


+ 18 


10 


10.0 


Ml 


140 


0.02 


3o,3 


-2o 


2, 4i,4 


-3i,3 




11(4) 


AA Tau'^ 


04 


34 


55.42 


+24 


28 


53.2 


K7 


140 


0.01 


3o,3 — 


2o,2 


. (4i,4 - 


-31,3)"^ 




<18 


LkCa 15" 


04 


39 


17.78 


+22 


21 


03.5 


K5 


140 


0.05 


3o,3 


-2o 


2, 4i,4 


-3l,3 




16(3) 


GM Aur<^ 


04 


55 


10.98 


+30 


21 


59.4 


K3 


140 


0.03 


3o,3 


-2„ 


2, 4i,4 


-3l,3 




12(5) 


IM Lup '^•'^ 


15 


56 


09.23 


-37 


56 


05.9 


MO 


150 


0.08 


3o,3 


-2o 


2, 4i,4 


-3i,3 




22(6) 


AS 209' 


16 


49 


15.29 


-14 


22 


08.6 


K5 


125 


0.028 


3o,3 — 


2o,2 


. (4i,4 - 


-31,3)"^ 




<30 


V4046 Sgrs 


18 


14 


10.47 


-32 


47 


34.5 


K5 


73 


0.01 


3l,2 


-2i 


1, 4i,4 


-3l,3 


~0.92j 


12(3) 


HD 142527'^ 


15 


56 


41.89 


-42 


19 


23.3 


F6 


145 


0.1 


3o,3 


-2o 


2, 4i.4 


-3l,3 




250(120) 



Star and dis k data from : lOi 



' Pinte et al 



'From this work and 



et al.l ll2004ll. lOi et al.l ll201ll). ^Andrews fc Williams! ll2005h. '^Lombardi et al.1 l|2008l 'l 



1 20081). 1 Andrews et al 



Oberg et al. 



1 201c 



I 2003), teodriguez et al. 



1 201 



Verhoeff et al. 



1 2OI1I) 



2011b ). Upper limits in parentheses. ^Calculated from other H2CO transition 
str engths and the determined excitation temperature. ''The 3 — 2 lines were reported as upper limits in 



Oberg et al. 



1 201c 



2011b), but after re-examination of the emission maps, we conclude that these lines are better described as tentative 

detections (2-3(t). 



